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We suggest a random field based model for calculation of physical properties of mixed ferroelectric 
relaxors. Our model naturally incorporates the different orientations of electric dipoles (related 
to different solid solution components) as well as the contribution of nonlinear and correlation 
effects of random field. We calculate the transition temperature Tc as well as concentrational 
and temperature dependence of order parameters. The equations for these quantities have been 
derived. The theory has been applied for description of the systems (PSN)i_a;(PST)a; with different 
degree of order as well as for (PMN)i_2,(PT)a; systems. We show, that higher Tc value for more 
disordered (PSN)i_i(PST)a; system at < a; < 0.5 is related to larger nonlinearity coefficient of 
PSN in comparison with that of PST. We determine the morphotropic region of temperatures and 
concentrations for (PMN)i_a;(PT)i,. The observed phase diagram of both aforementioned mixed 
ferroelectric relaxors is in pretty good coincidence with the results of our numerical calculations. 



I. INTRODUCTION 

The ferroelectric relaxors like PbMgi/3Nb2/3 03 
(PMN), PbSci/2Tai/203 (PST), PbSci/2Nbi/203 
(PSN), Pbi_^La:,Zro.65Tio.35 3 at a; = 0.08-0.09 (PLZT 
8-9/65/35) have been widely investigated in the last years 
[1]. In spite of this the mechanisms of unusual proper- 
ties of these materials (like Vogel-Fulcher (V-F) law in 
dynamic dielectric permittivity, nonergodic behaviour, 
the distribution of different properties maxima in broad 
temperature range (AT ^ 100 K for PMN relatively 
to the temperature of permittivity maximum)) are still 
under discussion. The presence of random field related 
to substitutional disorder, vacancies of lead and oxy- 
gen, impurity atoms is firmly established in the relaxors. 
Since this field has to influence the local properties of 
material, its distribution was shown to result in V-F law 
[2] , anomalies of nonlinear dielectric permittivity [3] , dis- 
tribution of relaxation times and non-Debye behaviour 
of dynamic dielectric permittivity [4] . The model [5] for 
quantitative description of relaxor ferroelectrics is based 
on the supposition that the random field destroys ferro- 
electric long range order which could appear at T < 
where is Burns temperature {T^ ^ 600 K for the most 
of relaxors). As a result, mixed ferroglass phase with 
coexistence of long and short range order or dipole glass 
state (both of them with non-ergodic behaviour) can ap- 
pear. In the case of PLZT the reference phase is PZT, 
its ferroelectric long range order is destroyed completely 



{a,t XLa — X > 0.07) or partly (at x < 0.07) by random 
field of impurity La ions. The degree of disorder can be 
controlled by special technological conditions which de- 
crease the number of random field sources in the relaxors 
PST, PSN and in other representatives of 1:1 family. In 
the case of PMN (1:2 family relaxor) it appeared im- 
possible to increase the degree of order with the help of 
technological treatment. However, this can be done by 
addition of PbTi03 (PT) (see e.g. [6] and ref. therein). 
In particular, measurements had shown that in solid 
solution (PMN)i_2;(PT)a; there is morphotropic phase 
transition between rhombohedral and tetragonal phases 
at Xc ~ 0.35. For x > Xc normal ferroelectric behaviour 
is present. Still missing theoretical description of this 
mixed system phase diagram seems to be necessary to 
clarify the physical nature and mechanisms of relaxor fer- 
roelectrics behaviour. On the other hand the solution of 
these problems is important due to giant electrostriction 
related to high dielectric response which result in ex- 
tremely high value of electromechanical coefficient (fc ^ 
0.92 for (PMN)o.9(PT)o.i single crystals). These anoma- 
lous properties made mixed ferroelectrics very attrac- 
tive for application in electromechanical transducers, 
actuators, sonars etc. [7,8]. Another mixed system 
(PSN)i_a;(PST)a; is also prospective for aforementioned 
applications [9]. Both phase diagram of these relaxors 
and the behaviour of PSN and PST with different degree 
of disorder exhibit the puzzles which have not been solved 
up to now. The most interesting problem is the increase 
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of transition temperature Tc in PSN and its decrease 
in PST as disorder increases. The observed dependence 
Tc{x) at < a; < 1 for solid solution of PSN and PST 
with different degree of order had shown [10,11,12] that 
aforementioned peculiarity exists in wide range of the 
components concentrations, namely T**(a;) > T°^'^{x) 
at X < 0.5. The experimental and theoretical investiga- 
tions of anomalies of dielectric response (the existence of 
low and high temperature maxima obeyed to V-F and 
Arrhenius law respectively) of PbSco.5(Nbo.2Tao.3)03 
single crystals have been performed recently [13]. In this 
work we propose the method of calculation of mixed sys- 
tems phase diagram. The method bases on the random 
field concept. We calculate the concentrational depen- 
dencies of transition temperature and order parameters. 
We show that the difference in PSN and PST behaviour 
is related to nonlinear and correlation effects contribu- 
tions of random field. The calculated phase diagrams fit 
pretty good the measured ones for (PMN)i_x(PT)j, and 
(PSN)i_,(PST),. 



II. THE MODEL 

Relaxors belong to the group of disordered ferro- 
electrics. In such materials random site electric dipoles 
tend to order system via indirect dipole-dipole interac- 
tion over soft mode of reference phase, while all other 
random electric field sources try to disorder the materials. 
The competition between these tendencies can lead to ap- 
pearance of coherently oriented dipoles. Their fraction L 
can be calculated in the random field theory framework 
on the base of equation 
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L= J (tjf{E,L)dE 



(1) 



where /(-B, L) is random field E distribution function. 
As the matter of fact L is a dimcnsionlcss order pa- 
rameter, which expressed through single dipole moment 



= {d* is effective dipole moment related to 

Lorenz field parameter, see [14] for the details) aver- 
aged over possible orientations and random field distri- 
bution function. In the case of two-orientation dipoles 
{h = IJcc = ly = 0) (Iz) = tanh{d* E;,/kT), Eq.(l) can 
be rewritten as 
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For the completely ordered system which is usually de- 
scribed by a mean field approximation f{E, L) = S{E — 
EqL), where Eq is mean field value {d*Eo = kTcmf), 
so that Eq.(2) gives well-known expression for order- 
disorder ferroelectrics: L^f = ta,nh{d* EoLmf /kT) and 

T —^ Tcmf at Lmf 0. 

For the considered case of mixed system of ferroelectric 
relaxors Ai-^^x two types of electric dipoles d^^ and 
and order parameters Li^ and Z/2z must be considered. 
In such a case the random field is induced by both type of 
dipoles and its distribution function F depends on both 
order parameters, i.e. F = f{E, Liz,L2z)- Therefore Liz 
and L2z can be written in the form of Eq.(2) with this 
distribution function and with dl^ and d2z in the Eqs. 
for Liz and L2z respectively, z and z being the directions 
of 1 and 2 types electric dipoles orientations. 

The most of the disordered system characteristics are 
included into random electric field distribution function 
f{E,L). This distribution function has been calculated 
in the statistical theory framework allowing for contri- 
bution of electric dipoles, point charges and dilatational 
centers as the random field sources [14]. Linear, non- 
linear and spatial correlation effects contribution of the 
random field was taken into account in [15]. Allowing for 
these effects for the materials with cubic symmetry (see 
[15]) in supposition that the considered electric dipoles 
are the main sources of random field and introducing 
cos(li^2) = cos(2:^;') one can obtain the following equa- 
tions for two order parameters: 
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Here a^''^^ is coefficient of nonlinearity of the reference phase, its dimensionality being s ^ (see [15]). 
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fi{ez,Li,L2) = ^ j exp{-plr]i)cos 
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f2{ez',Li,L2) = ^ j exp(-p^,972)cos 
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771 = C0S^(^1^2)Cl + 6; '72 = ^1 + 003^(^1^2)6; $1 
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In Eq.(Q) 6 and ^ describe the width of random field 
distribution functions with Gaussian form related, re- 
spectively, to the first and second types of dipoles, their 
concentrations being ni = and ri2 — '^^^V^'' (A is 
fraction of unit cells in which dipoles exist); £1,2, ai,2 and 
^■01,2 are respectively dielectric permittivity, lattice con- 
stants and correlation radii of A and B relaxors reference 
phase. Note, that the difference between distribution 
functions /i and /2 is related to the different orienta- 
tions of the dipoles which manifests itself in the presence 
of 003(^1^2) multipher; /i = /2 if cos{lil2) = 1- 

The polarization of mixed system P can be expressed 
via Li and L2 as follows 



(5) 



Integrals Ii and I2 depend on components concentra- 
tions since the distribution function width (see Eq. (||)) 
depends on x and Tc- 
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III. PHASE DIAGRAM OF MIXED SYSTEM 



The phase diagram has to describe the dependence of 
the temperature at which order parameters arise on 
concentration. The Tc{x) dependence can be obtained 
from Eqs.d), M in the limit Li ^ 0, L2 ^ 0. One 
can see from (jj), that in such a limit the distribution 
functions /i and /2 can be represented in the form: 



/i,2(ei, Li,L2) = /oi,2(£i) - fci,2 
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Substituting Eqs. (^)-(^) into (^ and performing the 
integration, one obtains the following system of equations 
for Tc{x) calculation: 
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FIG. 1. Concentrational dependence of dimensionless tran- 
sition temperature Tc = for the case di \\ d2, 
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a, 



(1) 




a, 



(2) 




1 (a): Q, 



(1) 




: a, 



(2) 




3(b). 



The condition of the system (H) solvability leads to 
following equation for Tc'. 



- C2{Tc)Te + Ci{Tc) = 0; 
Ci{Tc) = (1 - 2;)x/i/2r,„/iT„„/2(l - cos2(ir^)); 

C2{Tc) = xhTcrnfl + (1 " x)l2T, 



(11) 



cm/2- 
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The dependencies Ci^2(7c) are related to those of /1.2 
integrals (see Eqs.(||)), so that ( p^ is complicated non- 
linear equation for T^. It will be solved numerically. 
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Factor (1 — cos^(ZiZ2)) reflects the dependence on ori- 
entations of the dl and cia dipoles. It equals | in the case 
when the dipoles 1 and 2 are oriented, respectively, along 
[001] and [111] type of direction. This factor is zero as 
cos(Zi/2) = ±1 when the dipoles are parallel or antipar- 
allel to each other. In the latter case the system ( p^ ) 
reduces to single equation Tc = C2. One can see, that 
in general case there is at least two Tc values and the 
largest one has to describe the mixed system behaviour. 
It is seen from Eqs. (|ll|), ([l2|), (|), that Tc depends 
on several parameters - coefficients of nonlinearity, distri- 
bution function half- width and electric dipole moments 
ratio. To illustrate the influence of these parameters on 
Tc{x), we depicted in Fig.l the dimensionless transition 
temperature Tc = Tc{x) /Tcmf2 for the case when di \\ ^2 
(i.e. Tc = C2) for several values of dimensionless param- 
eters. It is seen that increase of aj, and A = d\/d2 
increases the rate of Tc{x) decrease. The influence of 
dipole moments ratio is more pronounced with decrease 
of the second component concentration. 
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FIG. 2. Phase diagram of mixed relaxors (PSN)i_^(PST):, 
for disordered (open symbols) and ordered (solid symbols) ce- 
ramic samples, taken from [10,11,12]. Solid and dashed lines 
- theory respectively for aforementioned and 10% larger A 
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disordered and ordered materials. 



IV. MIXED SYSTEM OF 1:1 FAMILY 
RELAXORS: (PSN)i x(PST)x 

The behaviour of this solid solution components (PSN 
and PST) is strongly different: the increase of the degree 
of disorder leads to increase of Tc for PSN and to de- 
crease of Tc for PST. At the first glance the latter seems 
to be more reasonable while the PSN behaviour is puz- 
zling since it is common wisdom that random field de- 
creases Tc- On the other hand the increase of random 
field value can result into appearance of nonlinear and 
correlation effects. These effects have to be dependent 



on the value of nonlinearity coefficient of the material. 
Keeping in mind that PSN and PST contain scandium, 
niobium and tantalum oxygen complexes, one can con- 
clude that coefficient of nonlinearity for PSN should be 
larger than that for PST. This is because this coefficient 
for niobium oxygen complexes was shown to be several 
times larger than that for tantalum [16]. The value of 
nonlinearity coefficients for PSN and PST was extracted 
from observed Tc values with the help of formulas (pd|), 
(|l|) at X = (PSN) and a; = 1 (PST). One can see, that 
in both cases Ci{Tc) = and 



Tc{PSN) = l2Tc^f2, 

Tc{PST) = hTc^fi. 



(13) 
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FIG. 3. Fraction of coherently ordered dipoles in 
(PSN)i-4PST),. 

a) sohd line - L2ord{x ~ 0); dashed line - L2disord{x = 0); dot- 
ted line - Liardix — 1); dotted-dashed hne - L\disord(x ~ 1) 

b) solid line - Pord{x — 0,3); dashed line - Pdiaord{x — 0,3); 
dotted line - Pord{x = 0,8); dotted-dashed line - 

Pdisord{x = 0, 8). 

The fitting of observed Tc values (see e.g. [10] and ref. 
therein) for more disordered PSN and PST made it possi- 
ble to obtain dimensionless coefficient of nonlinearity (ao 
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a3(fcT„„/2/rfI)) for PSN aj,' 



(2) 



0, 3 and a^^^ 



0,098 



for PST along with distribution function dimensionless 



half-width {qi^2 



_ 167r rft" /3 



C IDTT 



) 92 (PSN) 



= 0.425, 91 (PST) = 0.499. Note, that while introduc- 
ing dimensionless values we took into account, that pa- 
rameters of reference phase for PSN and PST are close 
to each other. Allowing for the fact that Tcm/ (PSN) 
~ 7cm/ (PST) [17] one comes to conclusion that ^ w 

d*(i2^ ~ which gives A^^ ~ 0.89 for more disordered 
relaxors. The same fitting of Tc values observed in more 
ordered PSN and PST [10,11] leads to = w 0, 
qi = 0.477, 52 = 0.416 and = 0.88. The ob- 

tained parameters made it possible to calculate Tc{x) 
both for more disordered and ordered (PSN)i_2;(PST)a; 
mixed system without additional fitting parameters. To 
simplify the calculations we supposed the same orienta- 
tions of di and d2 dipoles. The comparison of calculated 
and measured [10-12] Tc{x) values is shown in Fig. 2. We 
used two different values of A to show, that the concen- 
tration Xc at which Tc (disordered) — Tc (ordered) is sen- 
sitive to the parameter values. One can see that the the- 
ory describes the observed values pretty good, allowing 
for the distribution of measured Tc values related to dif- 
ferent experimental conditions (see [10] and ref. therein). 

The calculations of ii and L2 on the base of Eq.(||) 
had shown, that L2 > Li both for more disordered and 
ordered relaxors, the values of L for ordered ones be- 
ing larger than those for more disordered samples. The 
nonzero values of Li and L2 in whole range of the compo- 
nents concentrations give evidence that each component 
is in mixed ferro-glass phase (see also [13]). 



800r 

700- 
600- 
fc^ 500- 
400h 



300- 















cubic 




















yi 










1 

1 

1 


te 


tragon 


al 


pseudc 


1 

)-' 








cubic 


; 









0.2 



0.4 0.6 

X 



0.8 



FIG. 4. Phase diagram of (PMN)i_^(PT)^. Experimental 
data [6] - symbols; theory - solid line is Tc{x) and dashed line 
is morphotropic region. 

The results of numerical calculations for several con- 
centrations of solid solution components is depicted in 
Fig. 3. In Fig. 3a we represented the fraction of coher- 
ently ordered dipoles in PSN {x = 0) and PST {x = 1) 



which is proportional to these materials polarization. For 
two relaxors solid solution at a; = 0.3 and x = 0.8 we per- 
formed the calculations of dimensionless polarization (in 
the units of Po ~ ^ ) on the base of Eq. (^) (keeping in 
mind that di \\ c?2 and fii = [32) and showed the results 
in Fig. 3b. It is seen that increase of x leads to increase 
of polarization, P{x) of more ordered samples is larger 
than P{x) of disordered samples. 



V. MIXED SYSTEM (PMN)i x(PT)x 

This system is a solid solution of relaxor ferroelectric 
PMN (component 2) and normal ferroelectric PbTiOs 
(component 1). In fact, the phase diagram of this system 
can be calculated similar to that in the previous section, 
with some details different. 

First, due to the existence of normal ferroelectric phase 
transition in PbTiOs at — 763 K which can be cal- 
culated in a mean field approximation, Tcmfi — 763 K. 
Second, the distribution function of random field in or- 
dinary ferroelectrics is known to be (5-function because 
of negligibly small random field in the system, so that 
^1 = in Eq.(^). Note, that in such a case PT com- 
ponent contributes to random field distribution via PT 
mean field written as xTcmfiLi/d\ (see Eqs.(^). The 
number of electric dipoles ni = x/a^ where a is PT lat- 
tice constant. 
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FIG. 5. Temperature dependence of the fraction of coher- 
ently ordered dipoles in PMN (solid line) and PT (dashed 
line). 

It is important to emphasize that PMN relaxor state 
was considered as mixed ferroglass phase with coexis- 
tence of long and short range order. The peculiarities of 
nonlinear dielectric permittivity of PMN (see [18]) speak 
in favour of this statement. Keeping in mind that the 
dipoles in PMN and PT are oriented, respectively, along 
[111] and [100] type of directions, we calculated Tc{x) 
on the base of Eqs.(|l|), (H) with cosihh) = It 
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appeared possible to neglect the contributions of non- 
linear and correlation effects because they do not im- 
prove the fitting of calculated and measured Tc{x) depen- 
dence. The results of numerical calculations are reported 
in Fig. 4. Fitting parameters q2 = 0.51 and A = i were 
obtained from observed rc(PMN), Tc(PT) and = i, 

i.e. |i = so that only 11% of unit cells in PMN have 
electric dipoles because /3i — 1. One can see that our 
calculation describes pretty good observed Tc{x) depen- 
dence. 

To find the change of mixed system symmetry with 
variation of concentration of components, we performed 
the calculations of Li and L2 temperature and concentra- 
tional dependence. We start with calculations for PMN 
{x = 0) and PT {x = 1) using the aforementioned pa- 
rameters (see Fig. 5). The comparison of solid and dashed 
lines in Fig. 5 shows clearly the difference of the behaviour 
of order parameters in mixed ferroglass phase (PMN) and 
normal ferroelectric phase (PT) . The calculation of L2 (T) 
at T ^ had shown that L2(T = 0) « 0.58, so that the 
contribution of long range order is large enough in low 
temperature region in PMN. 



0.3 



,0.25 



0.2 



o 

•73 

N 

^0.15 



5 0.1 



























^ 


\ .\ 














■■.\\ 










A' 

\- 










.1 



^ 0.05[ 


0.2 0.4 0.6 0.8 1 
X 

FIG. 6. Temperature dependence of the components with 
rhombic (P2) and tetragonal (Pi) symmetry of dimensionless 
polarization in (PMN)i_^(PT):,. Solid line - P2{x = 0,1); 
dotted line - Pi{x = 0.3); dashed line - Pi{x = 0.1); dot- 
ted-dashed line - P2{x — 0.3). 

The polarization of mixed system was calculated with 
the help of Eq.(6), which incorporates the contribution of 
the first and second type of dipoles, i.e. P{x) = Pi{x) + 

P2{x). The concentrations at which Pi{x) = P2{x) 

correspond to coexistence of two types of symmetry in 
the mixed system, i.e. to appearance of morphotropic 
region in phase diagram. Temperature dependence of di- 



Pi 



and 



P2 



for several x values is depicted 



mensionless 
in Fig. 6. 

The most interesting feature in the Fig. 6 is the ex- 
istence of the curves crossover at some T = Tcr{x) for 



X = 0.3. This means that |Pi| = jPal at T = Tcr{x) 
(morphotropic region), |Pi| > jPj] at T > Tcr{x) or 
|Pi| < IP2I at T < Tcr{x), i.e. mixed system has the 
symmetry of component 1 or component 2 respectively. 
For (PMN)i_^(PT)^ the value T^rix = 0.3) = 275 K and 
Tcr{x — 0.4) — 460 K, therefore Tcr increase with concen- 
tration increase (see dashed line in Fig. 4). On the base 
of this consideration we discern the symmetry of mixed 
(PMN)i_a;(PT)j, system in Fig.4 (pseudo-cubic like PMN 
and tetragonal like PT) which is in agreement with ob- 
served one [6]. 

One can see from Fig. 6 the absence of the curves 
crossover, i.e. morphotropic region at x = 0.1. Since 
the compositions with small concentration of PT (x « 
0.1) are important for applications (see [7,8]) we de- 
picted in Fig. 7 the fraction of coherently ordered dipoles 
and the components of polarization with tetragonal 
(solid line) and pseudo-cubic (rhombohedral) symmetry 
(dashed fine) at T = 0. 
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FIG. 7. Concentrational dependence of order parameters 
at T = with pseudo-cubic (rhombohedral) (dashed lines) 
and tetragonal (solid lines) symmetry expressed via fraction 
of coherently ordered dipoles (a) and via contributions to po- 
larization (b). 
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It is seen from Fig. 7b that tetragonal component lin- 
early increases with PT concentration increase whereas 
pseudo-cubic component has a maximum at a; w 0.135. 
This maximum origin can be related to the competition 
between increase of L2 (see dashed line in Fig. 7a) and 
(1 — x) decrease with x increase (see Eq.(6)). Because of 
pseudo-cubic symmetry, both the polarization maximum 
and maxima of dielectric response, piezoelectric and elec- 
tromechanical coefficient can be expected at x = 0.135. 
This concentration is a little larger than x w 0.1 where 
the high values of aforementioned properties were ob- 
served [7,8]. 



VI. CONCLUSION 

We propose a model for calculation of ferroelectric or- 
der parameters and phase diagram of mixed relaxors. 
The physical background of the model is influence of 
random electric field of mixed system on its properties. 
Randomly distributed electric dipoles were supposed to 
be the main sources of the field. The contribution of non- 
linear and spatial correlation e£Fec;ts of random field was 
taken into account. We carried out the specific calcula- 
tions for mixed PSN-PST and PMN-PT systems. The 
solution of the puzzle of larger transition temperature 
of more disordered (PSN)i_j:(PST)a; system in the re- 
gion < X < 0,5 was shown to be related to nonlin- 
ear and correlation effects contribution which has to be 
larger in PSN. The obtained concentrational dependence 
of transition temperature and polarization of solid so- 
lution (PMN)i_j;(PT)j; revealed the existence of mor- 
photropic region with coexistence of rhombic (pseudo- 
cubic) and tetragonal symmetry phases, its concentration 
being dependent on temperature. The maximal contribu- 
tion of pseudo-cubic symmetry polarization at x = 0.135 
was supposed to be related to high values of electrome- 
chanical coupling coefficient, observed at x w 0.1 in the 
mixed system. The developed theory describes pretty 



good observed phase diagram of PSN-PST and PMN-PT 
mixed relaxor systems. 
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